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ABSTRACT 

The Survey of Hi in Extremely Low-mass Dwarfs (SHIELD) is an on-going multi-wavelength pro¬ 
gram to characterize the gas, star formation, and evolution in gas-rich, very low-mass galaxies that 
populate the faint end of the galaxy luminosity function. The galaxies were selected from the first 
~ 10% of the Hi ALFALFA survey based on their low Hi mass and low baryonic mass. Here, we 
measure the star-formation properties from optically resolved stellar populations for 12 galaxies using 
a color-magnitude diagram fitting technique. We derive lifetime average star-formation rates (SFRs), 
recent SFRs, stellar masses, and gas fractions. Overall, the recent SFRs are comparable to the life¬ 
time SFRs with mean birthrate parameter of 1.4, with a surprisingly narrow standard deviation of 
0.7. Two galaxies are classified as dwarf transition galaxies (dTrans). These dTrans systems have 
star-formation and gas properties consistent with the rest of the sample, in agreement with previous 
results that some dTrans galaxies may simply be low-luminosity dlrrs. We do not find a correla¬ 
tion between the recent star-formation activity and the distance to the nearest neighboring galaxy, 
suggesting that the star-formation process is not driven by gravitational interactions, but regulated 
internally. Further, we find a broadening in the star-formation and gas properties (i.e., specific SFRs, 
stellar masses, and gas fractions) compared to the generally tight correlation found in more massive 
galaxies. Overall, the star-formation and gas properties indicate these very low-mass galaxies host a 
fluctuating, non-deterministic, and inefficient star-formation process. 

Subject headings: galaxies: dwarf - galaxies: evolution - stars: Hertzsprung-Russell diagram 


1. INTRODUCTION 

The evolution of a galaxy is often strongly affected 
by the star-formation process. Numerous studies have 
characterized the relationship between star-formation ac¬ 
tivity, gas conditi ons, local environment, and galaxy 
morpholog y (e.g.. iKennicuttl 1199ft Balogh et all 12004 
Kennicutt fc Evandl2012l).~Vohime- limited surveys (e.g., 
Brinchmann et al. 12004 iLee et all 120071 : iBothwell et al'l 
20091) have mapped galaxy star-formation properties over 
stellar mass scales of 7 < log (M*/M 0 ) <11. These sur¬ 
veys have found that the properties vary monotonically 
over a luminosity range —19 <Mg < —15 mag, but the 
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distribution of star forming properties widens at fainter 
magnitudes roughly corresponding to a stellar mass of 
log (M*/M q ) of 8.5. Specifically, galaxies at the lower 
mass ranges of these larger surveys show a broadening in 
their star-formation and gas properties, with less correla¬ 
tion between gas mass, stellar mass, and star-formation 
activity. 

While these large surveys, as well as smaller, more 
detailed studies (e.g., VLA-ANGST; lOtt et all 12012 ). 
have begun characterizing the star-formation and gas 
properties (and their dispersion) in lower mass galaxies, 
the faint end of the galaxy luminosity function remains 
largely unexplored. Outside of a dense group or clus¬ 
ter environment, these dwarf galaxies are expected to 
be gas rich, presumably evolving in relative quiescence 
due to inefficient star formation. Generally thought to 
have longer gas-cooling timescales and low, if not sup¬ 
pressed, star-formation activity, they play a key role in 
disentangling the factors driving star formation and star- 
formation driven galaxy evolution. Further highlight¬ 
ing the importance of low-mass galaxies, recent studies 
show that low-mass galaxies contribute significantly to 
both the total sta r-formation rate density at high redshift 
(e.g. JAla vi et al.||2014D and the ionizing backgro und ra¬ 
diation (jMostardi et al.ll2013t iNestor et al.ll2013H . While 
galaxies at the faint end of the luminosity function are 
predicted to be the most numerous class of galaxies, their 
intrinsic low luminosities and small angular sizes have 
made identifying and cataloging these systems through 
optical surveys difficult. Their high gas fractions make 
them better candidates for discovery through Hi surveys. 

The Arecibo Legacy Fast ALFA (ALFALFA) survey 
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(|Giovanelli et al Jp2005l : lHavnes et al.ll20lll) is a blind ex- 
tragalactic Hi survey to map the nearby Hi universe over 
7000 deg 2 of high Galactic latitude sky. With an Hi mass 
detection limit as low as 10 6 M 0 for galaxies in the Lo¬ 
cal Group and 10 9 5 M 0 at the survey velocity limit of 
2 ~ 0.06, the ALFALFA survey was designed to populate 
the faint end of the Hi mass function over a cosmologi- 
cally significant volume. 

The Sur vey of Hi in Extre mely Low mass Dwarfs 
(SHIELD; ' Camion etaD l20ll was designed to exam¬ 
ine the early release of the ALFALFA dataset in a sys¬ 
tematic investigation of nearby galaxies with HI masses 
< 10 7 M 0 outside the Local Group. From the first 
~ 10% of the processed ALFALFA survey data, 12 
systems were selected for further study. HST opti¬ 
cal imaging of these 12 systems have facilitated accu¬ 
rate distances bas ed on the tip of the r ed giant branch 
(TRGB) method (jMcQuinn et al.l 120 1 11 1 . The distances 
range from 5 to 12 Mpc; the Hi masses range from 
4 x 10 6 - 6 x 10 7 M 0 based on these distances. All 
of the CMDs have red giant branch (RGB) sequences; 
a subset of the systems also show evidence of recent 
star formation based on their more populated upper 
main sequences (MS) and helium burning se quences 
(|Dohm-Palmer et al.l 119971 : iMcQuinn et all 1201 If) . Most 
of the galaxies are detected in Ha and spectroscopy of 
their Hll re gions reveal that t h e SHIEL D galaxies are 
metal poor (lHaurberg. Salzer. fc Cannonll201 4). 

In addition to providing accurate distance measure¬ 
ments, resolved stellar populations provide a means to 
measure the star-forming properties of a system by 
reconstructing a star formation history (SFH) using 


gram fCMD) fitting 

techniques (e.g., Tosi et al. 

1989 

Tolstov & Saha 

1996; 

Gallart et al. 19961: Dolphin 

1997; 

Holtzman et al. 

1999; 

Harris & Zaritskv 200lf). This ap- 


proach has been used to measure the star-formation rate 
as a function of time (SFR(t)) in a significant num¬ 
ber of nearby galaxies, putting constraints on star for- 
mati on, star-formation feedback, and galaxy evolution 
fe.gJGa llart fo LCID Teainll2007t iDalcanton et, al.l[2009f 

McQ uinn et al.ll2010afl . 


Here, we use the optically resolved imaging of the 
SHIELD galaxies to measure the star-formation charac¬ 
teristics of gas-rich galaxies in a mass range that is not 
well probed in previous studies. In this work, we use 
t he photometry from HST optical images measured in 
IMcQuinn et al.1 ( 2014 ) to derive the star-formation prop¬ 
erties using a CMD fitting technique and stellar evolution 
isochrones. In Section 2, we describe the sample, obser¬ 
vations, and data reduction. In Section 3, we explain 
the CMD fitting technique used to measure the star- 
formation activity. In Sections 4 and 5, we present the 
lifetime average SFRs, recent SFRs, and discuss whether 
the star-formation activity appears to be correlated with 
the local environment. In Section 6, we compare the 
star-formation and stellar mass properties with the_Hi 
gas pr operties previously measured from iCannon et all 
( 2011 ?) and compare the results in these low-mass galax¬ 
ies to previously published works. Finally, in Section 7 
we summarize our conclusions. 


2. THE GALAXY SAMPLE, OBSERVATIONS, AND 
PHOTOMETRY METHOD 


~i i i i i i r 



Figure 1 . HST composite image of AGC 110482, representative 
of the SHIELD sample, combining the F606W image (red), the av¬ 
erage of F606W and F814W images (green), and the F814W image 
(blue). North is up and east is left. The field of view encompasses 
twice the optical major axis as determined by iteratively plotting 
CMDs with different elliptical parameters. 


Table □ lists the 12 galaxies that comprise the 
SHIELD sample, their coordinates, observation details, 
and TRGB distances. The observations consist of optical 
imaging from the Adv anced Ca m era f or Surveys (ACS) 
Wide Field Channel dFord et al.lI1998T) aboard the HST 
in the F606W and F814W filters. New observations of 
11 galaxies were obtained as part of the HST-GO-12658 
SHIELD program (PI: Cannon). Archival observations 
were available for three galaxies from the previous HST- 
GO-10210 program (PI: Tully); two of these targets over¬ 
lap with our new observations providing additional pho¬ 
tometric depth to the final images. All observations were 
obtained with two filters during single HST orbits and 
were cosmic-ray split (CRSPLIT=2). Average integra¬ 
tion times were 1000 s in the F606W filter and 1200 s in 
the F814W filter. The observations for the two galaxies 
that overlapped in the observing programs have longer fi¬ 
nal integration times of ~1900 s in the F606W filter and 
~2300 s in the F814W filter. The observing programs 
and final integration times for each galaxy are noted in 
Table □ 

Figure [T] shows an example composite image from the 
data of AGC 110482, reproduced from the previous st udy 
meas uring the distances to the galaxies (IMcQuinn et al.1 
12014 see their Figure 1). The images were cosmic-ray 
cleaned, processed by the standard ACS pipeline, cor¬ 
rected for charge transfer efficiency non-linearities, and 
combined using MultiDrizzle. The F606W filter image 
is shown in blue, F814W in red, and the average of the 
F606W and F814W images in green. Typical of the sam¬ 
ple, and of gas-rich dwarf galaxies in general, the image 
shows both a young, blue stellar population and a red, 
older and somewhat more extended stellar population. 

The sample has also been imaged in the narrow- 
band Ha filter from th e WIYN 3.5 m telescope 
(lHaurberg. Salzer. fe Cannonl 12014 . These images re¬ 
veal that 2 of the 12 galaxies lack Hu regions. Specif¬ 
ically, AGC 749241 was not detected in Ha and 
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Table 1 

Summary of the SHIELD Sample and Observations 




R.A. 

Decl. 

HST 

F606W 

F814W 

TRGB Dist 

+606W 

Galaxy 

Alt. Name 

J2000 

J2000 

ID 

(sec) 

(sec) 

(Mpc) 

(mag) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

AGC 110482 

KK 13 

01:42:17.4 

26:22:00 

12658 

1008 

1226 

7.82i0.21 

0.26 

AGC 111164 

KK 17 

02:00:10.1 

28:49:52 

10210 

934 

1226 

5.11T0.07 

0.16 

AGC 111946 

KK 15 

01:46:42.2 

26:48:05 

12658 

1008 

1126 

9 02” 1 ” 0 ' 20 
-0.29 

0.24 

AGC 111977 

KK 16 

01:55:20.2 

27:57:14 

12658 

1008 

1226 

5 96 -1 ” 0 ' 11 
o.yo_o 09 

0.20 





10210 

934 

1226 


AGC 112521 


01:41:07.6 

27:19:24 

12658 

1008 

1226 

6.58T0.18 

0.18 

AGC 174585 


07:36:10.3 

09:59:11 

12658 

1000 

1220 

7.89i°;?i 

0.11 





10210 

900 

900 


AGC 174605 


07:50:21.7 

07:47:40 

12658 

1000 

1200 

10.89i0.28 

0.07 

AGC 182595 


08:51:12.1 

27:52:48 

12658 

1008 

1226 

9.02i0.28 

0.12 

AGC 731457 


10:31:55.8 

28:01:33 

12658 

1008 

1226 

11 10+0.20 
11.1«3_0 16 

0.08 

AGC 748778 


00:06:34.3 

15:30:39 

12658 

1008 

1226 

a Aa + 0.14 

6.46I 0 it 

0.19 

AGC 749237 


12:26:23.4 

27:44:44 

12658 

1008 

1226 

11 62~l” 0 ' 20 

11 .oz—o 16 

0.06 

AGC 749241 


12:40:01.7 

26:19:19 

12658 

1008 

1226 

5 62”*~ 0 ' 17 

°- DZ —0.14 

0.04 


Note. — Column 1—Galaxy name. Column 2—Alternate galaxy name. Columns 3 and 4—Coordinates 
of galaxy in J2000. Column 5—HST observing program. Columns 6 and 7—Int egration time in the F 606W 
and F814W filters with the ACS instrument. Column 8—TRGB Distances from [McQuinn et al ] mm- The 
distances were derived using a Maximum-Likelihood technique; the uncertainties re flect the asymmetric d istri- 
bution of photometric errors. Column 9—Galactic a bsorption from the dust maps of ISchlegel et aL| d 1998T) with 
the recalibration from Schlafly & Finkbeiner (2011). 

AGC 748778 had only weak, diffuse Ha emission but 
no Hii region. These Ha data suggest that the star- 
formation activity over the most recent ~ 5 Myr has 
been either below detection or very low (Ha-based SFR 
< HT 5 M 0 yr” 1 ). Thus, AGC 749241 and AGC 748778 
meet the criteria of having both detectable amounts of Hi 
and no significant recent SFR as measured by Ha and can 
be classified as transition dwarfs (dTr ans; iMateol 119981 
iSkillman et al J 120031 : 1 Weisz et alJl2011[ ). Ha-based SFRs 
for the rest of the sample will be presented in N. C. Hau- 

rberg (in preparation). _ _ 

As described in detail in [McQuinn et al.l (120141) . 
photometry was performed on the pipeline processed, 
cosmic-ray cleaned images (CRJ files) using the ACS 
module of the DOLPHOT photometry package (|Dolphinl 
1200(f ). Both quality and spatial cuts were performed on 
the photometry. Briefly, the photometry was filtered for 
well-fit point sources recovered with a signal-to-noise ra¬ 
tio > 5 that were not severely affected by stellar crowd¬ 
ing. Spatial cuts were determined iteratively by plotting 
the CMD of stars in concentric ellipses with the same 
ellipticity and position angles centered on the optical 
center of each galaxy. In the inner regions, the CMDs 
are dominated by stars from the galaxies, while in the 
outer regions, the CMDs are dominated by background 
sources. The semi-major and semi-minor axes of the el¬ 
lipses were increased until the CMDs from larger annuli 
matched the distribution of point sources from a field re¬ 
gion CMD. Artificial star tests using ~500k stars were 
performed to measure the completeness limit of the im¬ 
ages using the same photometry package and filtered on 
the same parameters. 

Figure [5] presents an example CMD for one of the far¬ 
thest galaxies from the SHIELD sample, AGC 749237, 
plotted to the 50% completeness level as determined by 
the artificial star tests, and one of the closest galaxies, 

AGC 111977. Representative photometric uncertainties 
per magnitude are shown for both. The photometry in 


the CMDs were cor rected for Galactic a bsorption based 
on the du st maps o f iSclilegel et all d!998l) with recalibra¬ 
tion from iSchlaflv fc Finkbeined ( 20111) : these values are 
also noted in Table |T| The depth of photometry shown 
in Figure [2] brackets that of the overall sample. For the 5 
galaxies that are located between 8 — 12 Mpc, the depth 
of photometry is similar to AGC 749237, reaching ap¬ 
proximately 1 mag below the TRGB, but populated by 
differing numbers of stars. For the 7 galaxies that are 
located between 5 — 8 Mpc, the depth of photometry 
is similar to AGC 111977, reaching approximately 2 mag 
below the TRGB. CM Ds of all 12 galaxies can be found in 
iMcQuinn et al.l (|20l4 see their Figure 2). In each CMD, 
the MS is identifiable as well as a populated RGB se¬ 
quence. The presence of stars in the upper MS indicates 
that a system has experienced star formation at recent 
times (t< few 100 Myr), while the presence of the RGB 
stars indicates star-formation activity at older times 
(t> 1 Gyr). Some of CMDs also show somewhat popu¬ 
lated helium burning sequences, including AGC 110482, 
AGC 111977, AGC 731457, AGC 749237, which is a fur¬ 
ther indication of recent star-formation activity. Asymp¬ 
totic giant branch (AGB) populations, indicative of in¬ 
termediate age star-formation activity, are also observed 
in a subset of the CMDs, most notably in AGC 110482, 
AGC 174605, and AGC 182595. 


3. METHODOLOGY FOR DERIVING THE 
STAR-FORMATION PROPERTIES 

The star-formation characteristics of our sample 
of galaxies are reconstructed using stellar evolution 
isochrones to create a series of modeled CMDs from syn¬ 
thetic stellar populations of different ages and metal- 
licities. The SFH of the best-fit modeled CMD to the 
observed CMD represents the most likely star-forming 
properties of the system. We use the numerical CMD 
fitting program, MATCH (iDolphinl 120021) . w ith the stel¬ 
lar evolutionary models from iMarigo et aD (I2008D and 
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Figure 2. CMDs of AGC 749237 (left) and AGC 111977 (right). 
Average uncertainties per magnitude are shown in each panel. The 
photometric depth of the CMD for AGC 749237 is typical of the 5 
farthest galaxies, while the depth of the CMD for AGC 111977 is 
typical of the 7 closest galaxies. The number of stars in the CMDs 
across the sample varies between ~ 800 and ~ 6500. 

updated AGB tracks from lGirardi et ahl (12010Q . The pri¬ 
mary inputs in reconstructing the past SFRs are the pho¬ 
tometry and artificial star recovery fractions (i.e., quan¬ 
titative observational uncertainties and incompleteness) 
coupled with the stellar evolutionary models. A Salpeter 
single sloped power law initial mass fun ction with a spec- 
tral index of —1.35 from 0.1 — 120 M 0 ; Salpete r 10-55 ) is 
assumed as well as a binary fraction of 35% with a flat 
secondary mass distribution. Distance is a free parame¬ 
ter fit by the SFH recovery program. We constrained 
the mean metallicity, Z(t), to be a continuous, non¬ 
decreasing function with time. This physically motivated 
constraint guides the metallicity evolution in the absence 
of observational constraints that would be available with 
deeper photometric data and produces a more realistic 
metallicity evolution in the galaxies avoiding large jumps 
in metallicity over short time periods. 

Photometric errors and extinction can broaden fea¬ 
tures in a CMD and are explicitly accounted for in the 
CMD fitting program. The photometric uncertainties 
and completeness are quantified using the artificial star 
tests. Extinction is a free parameter fit by the SFH 
recovery program. Thus, while the photometry shown 
in the CMDs in Figure [2] has been corrected for fore¬ 
ground extinction, we used the uncorrected photometry 
in the CMD fitting program. Both foreground and in¬ 
ternal extinction are expected to be low for this sample. 
The galaxies are located at high Galactic latitudes with 
foreground extinction estimated to range from Apeo6W 
of 0.04 to 0.26 mag (see Table [TJ. Internal extinction 
is also expected to be low as dwarf galaxies ha ve been 
show n to follow the mass-metallicity relation dBerg et al.l 
HMD. Nonetheless, similar to the distance, extinction is 
an important parameter to fit in measuring the SFHs 
and provides a consistency check when compared to in¬ 
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dependently measured values. For our sample the dis¬ 
tances and extinctions derived from the best-fit CMDs 
were within 0.1 mag and 0.05 mag, respectively, of the 
measured TRGB distances and foreground extinction es¬ 
timates. Once the best fit has been verified, the final 
SFH solution is derived with the distance and extinction 
parameters fixed to the measured values from Table [T] 

Figure [3] presents two examples of the observed CMDs 
alongside the modeled CMD on the same axis scale for 
comparison as well as the residuals and weighted resid¬ 
uals between the observed and modeled CMD. As seen 
in Figure [U the models characterize all of the features 
of the CMDs quite well with only small discrepancies. 
The synthetic CMD of the galaxy is based on the most 
likely SFH and metalli city of the ga laxy g iven our in¬ 
puts and models (e.g.. iDolphin et aLl 12003T ). While the 
constraints on the metallicity evolution of the sample are 
limited by the photometric depth of the data, the best-fit 
values of the most-recent metallicities are in agreement 
with spectroscopic oxygen abundance measurements and 
the mass-metallicit y relati on for l ow-ma ss galaxies as dis¬ 
cussed in [Haurberg. Salzer. fe Cannonl d2014l) , providing 
an additional consistency check on our results. 

Uncertainties on the SFHs take into account systematic 
and random uncertainties. The systematic uncertainties 
measure the uncertainty in the SFHs due to variations of 
the real data fr om the theoretical stellar evolution models 
(jDolphinil2012l) . These uncertainties are estimated by ap¬ 
plying shifts to the models in both bolometric luminosity 
and temperature using Monte Carlo simulations. Ran¬ 
dom uncertainties were estimated using a new method of 
applying a hybr id Marko v Chain Monte Carlo simulation 
as described bv iDolphlnl (12013D . This technique samples 
the probability distribution of the parameters used in the 
original solution to directly estimate confidence intervals. 
The inclusion of this additional source of uncertainty is 
particularly important for time bins where the SFR is 
low or zero, and for data whose photometric depth is 
above the horizontal branch or red clump as is the case 
for the current galaxy sample. The final uncertainties in 
the SFHs are calculated by adding the uncertainties in 
quadrature. 

Reconstructing detailed SFHs (i.e., SFRs(t)) requires 
sufficient information on the ages and metallicity of the 
stellar populations from a CMD. Therefore, the photo¬ 
metric depth of the data is the main determinant on the 
temporal resolution achievable in a SFH at older ages. 
For a discussion on the systema tic effects of pho t ometri c 
depth on measuring a SFH see iMcQuinn et ahl (|2010aD . 
A secondary determinant on the temporal resolution is 
the number of stars populating a CMD, which becomes 
a more significant factor to consider in low-mass galax¬ 
ies. The small number of stars in the SHIELD galaxies 
compare d to m ore mas sive dwarfs (jDalcanton et al.ll2009l : 
IMcQuinn et al.ll2010al) is therefore one of the dominant 
limitations in determining recent SFHs with high time 
resolution. Through extensive tests of MATCH using 
different time binning schemes, we found that measur¬ 
ing the SFRs over two timescales, recent (t~200 Myr; 
primarily measured by the upper part of an optical 
CMD) and lifetime-average SFRs (primarily measured 
from the RGB stars), resulted in the most robustly mea¬ 
sured SFRs across the sample. The recent timescale of 
200 Myr is an appropriate choice for the entire sample, 
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(b) modeled CMD 


(a) observed CMD 


(b) modeled CMD 



Figure 3. Left 4 panels: AGC 749237: Hess diagrams of the observed CMD compared with the best-fit modeled CMD, residual of the 
data—model CMD with bl ack and white p oints representing iL5cr, and the residual significance CMD which weights the data—model by the 
variance in each Hess bin dDolphinl 12002T) . The synthetic CMD reconstructs the different features seen in the observed CMD with minor 
discrepancies. The axes labels combine the name of the HST instrument (ACS WFC) with filter number. Right 4 panels: AGC 111977: 
Same sequence of hess diagrams as shown for AGC 749237. 


regardless of distance and photometric depth, because 
of the generally low star-formation activity in low-mass 
galaxies. Further, the ~200 Myr timescale for the recent 
SFRs is comparable to the timescale measured by near 
UV studies of star-for mation rates^ in particular from 
GALEX images (e.g., IHao et all 1201 111 . Thus, the re¬ 
cent SFR measurements reported here will enable com¬ 
parisons with a large number of galaxies that have SFRs 
measured from integrated UV light. 

Note that the data for a few of the galaxies allow only 
upper limits to be placed on the SFRs. For the lifetime 
SFRs, this is mainly due to limits in photometric depth 
and impacts two of the three farthest galaxies. For the 
more recent SFRs, this is a result of the low number 
of young stars available to fit in the CMD and impacts 
the two galaxies with the lowest recent SFRs. The re¬ 
cent and lifetime SFRs, as well as the uncertainties, are 
reported in Table [2] For the galaxies whose measured 
SFRs are upper limits, we report only the upper bound 
of the uncertainties. 

Stellar mass estimates can be derived from the lifetime- 
average SFRs. Integrating the lifetime-average SFR over 
the life of a galaxy yields the total stellar mass created 
in a system. This total stellar mass can be converted to 
a measurement of the current stellar mass by applying a 
correction for the amount of mass returned to a galaxy 
over the lifetime of th e stellar populati o ns. W e used the 
recycling model from iKennicutt et all (11994T) who esti¬ 
mate that 30% of the mass in a stellar generation is re¬ 
turned to a galaxy over its lifetime based on a Salpeter 
IMF. These values are reported in Table [5] As a consis¬ 
tency check on our stellar mass values, we compared our 
results with preliminary estimates of the stellar masses 
based on mass-to-light (M/L) ratio analysis using Spitzer 


Space Telescope IRAC imaging (J. M. Cannon et al., in 
preparation). Our stellar mass estimates are a factor of 
2 higher than the majority of the stellar mass estimated 
from the M/L ratios, but are still in general agreement 
with one another given the large uncertainties. The ex¬ 
ception is AGC 748778 whose stellar mass estimate from 
the M/L ratio is significantly lower than our value. Fu¬ 
ture work on the SHIELD sample will include a larger 
discussion of the two stellar mass estimates including a 
description of the infrared imaging, analysis, and uncer¬ 
tainties (J. M. Cannon et al., in preparation). 

4. STAR-FORMATION RATES AND COMPARISON WITH 
OTHER SURVEYS 

Figure [I] presents a comparison between the recent 
(t< 200 Myr) SFRs with the average lifetime SFRs for 
the SHIELD sample, dlrrs are plotted as filled, blue cirl- 
ces; dTrans are plotted as unfilled, black cirlces. The 
recent and lifetime SFRs are low, ranging from —3.5 < 
log (SFR) < —1.8, with SFRs in units of M 0 yr -1 . The 
top axis in Figured shows the stellar mass for each sys¬ 
tem, which provides a complementary view of the re¬ 
sults. Seen in this way, it is clear that the galaxies follow 
the well-established correlation that higher mass galaxies 
have higher recent SFRs. This trend reflects that while 
gas fractions tend to decline with increasing galaxy mass, 
total gas masses increase as a function of mass and can 
thus support higher SFRs (e.g.. iBoselli et al.1l200lH . 

In Figure [2 the solid line marks where the recent to 
lifetime SFR ratio (i.e., the birthrate parameter b = 
SFR recent / SFR Hfetime ; [Scaloj[i986D is equal to 1. The 
birthrate parameters measured for the sample are listed 
in Table [2] The range in the b parameter is narrow with 
a mean value of 1.4 and a dispersion of 0.7. To test 
whether the dispersion is due to measurement errors or 
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Table 2 

Star Formation and Gas Properties 



SFR 

SFR 

Recent / 

Total 


Gas 



< 200 Myr > 

< Lifetime > 

Lifetime 

M* 

M HI 

Ratio 

T gas 

Galaxy 

(10- 3 M 0 yr- 1 ) 

(10- 3 M 0 yr- 1 ) 

SFR 

(10® M 0 ) 

(10® M 0 ) 

(Mhj/M.) 

(Gyr) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

AGC 110482 

5 4 +2-3 
-2.2 

5.5 + 1 '® 

1.0 ±0.5 

55 ±\ 9 g 

19.±2 

0.35±0.13 

6.7±2.9 

AGC 111164 

0 52 ~*~ 0,23 
u oz -0.35 

1 Q +0-2 
i U -0.3 

Q c +0.3 
U -° -0.4 

10 ± 2 

4.0±0.4 

0 42 - '" 0,09 
u -^ z -0.15 

14+ 7 

x ^-10 

AGC 111946 

o c +1-7 

0.0 —2.1 

i n +0.6 

-0.7 

2 1 + 1,3 
z± -1.5 

17 t 6 7 

15±1.7 

n ss+0-35 

U.88-0 38 

17 q+4.0 
7.o_4 ? 

AGC 111977 

2.6 j),- 3 

3.8 ±l\l 

0 7 + 04 
u -' -0.4 

07 +12 

01 -11 

7.1±0.8 

0.19±0.06 

r 0 + 2.7 

0.0_2 4 

AGC 112521 

< 0.36+ 014 

0 75 + 0 - 26 

U -'° -0.25 

< 0.5+ 0 3 

7 +3 

1 -2 

7.0±0.8 

0.94j)°; 3 ® 

> 37+ 15 

AGC 174585 

n o +1.2 

-1.3 

0 90 + 0,33 
u.uu _ 0 30 

2 5 + 1 ' 6 
-1.6 

9 ”*” 3 
y -3 

7.9±0.9 

0.89l°; 34 

6-6±l;? 

AGC 174605 

4.3 + \l 

< 2.8+ 14 

> 1.5+ 10 

< 28+ 14 

19±2 

> 0.66+ 0 33 

8.2 +i A g 

AGC 182595 

4.8 i 3 ;° 

r i +2.3 

0,1 -3.3 

0 Q +0-7 
uy -0.6 

50 tf 2 

8.1±0.9 

o.i6t°;?o 

O 9 + 2.1 

°’ z — 0.9 

AGC 731457 

14. ±6. 

6 6 ”*~ 3 ' 7 
o.o _ 4 g 

2 1 +1,5 
z± -1.8 

65 til 

18±2 

0.281®; 1 ® 

2.5±1.1 

AGC 748778 

0-68 t° 0 fs 

q 09 +0.06 
u-oz _ 0 -q 

9 -1 +0.9 
-1.4 

3 

4.5±0.5 

(M 

CO IO 

6 d 
+ 1 

13±® 

AGC 749237 

i3. +;■ 

< 5.4+ 2 - 9 

> 2.4+ 1 - 8 

< 53+ 29 

57±6 

> 1 ^+0.60 

8.4+ 4 -® 

AGC 749241 

< 0.32+ 0 08 

0 36 + 0 13 

U.OO _0.24 

< 0.9+ 0 4 

4 " l ” 1 
^ -2 

5.7±0.7 

1 f?+0.60 
± -°-l.l 

> 34+ 9 


Note. — Column 1—Galaxy name. Columns 2—3 Average SFRs over the past 200 Myr and the lifetime of the galaxy 
derived from the CMD fitting technique. Uncertainties include both systematic and random uncertainties; see text for 
details. For cases where only an upper limit on the SFR could be determined, no lower bound on the uncertainty is listed. 
Column 4—Ratio of the recent to lifetime average SFR. Columns 5— Stellar Masses formed over the lifetime of the galaxies, 
calculated from the SF Rs in column 3 and assuming 30% of the stellar mass is re cycled over th e lifeti me of the galaxy 
JKennic utt ct al. 1994) . Column 6—Hi mass based on the Hi flux measurements from [Cannon et al 3 1 2013) and the adopted 
TRGB distances from iMcQuinn et af ] (j2014b . Note the Hi masses reported for AGC 112521 and AGC 182595 have been 
updated. Columns 7— Gas Ratio (M#/ / M*). Column 8—Gas consumption timescales based on the recent SFRs, assuming 
a 30% recycling fraction of mass, and the Hi masses with a 1.33 correction for helium. 


log (M./M 0 ) 

6.0 6.5 7.0 7.5 8.0 8.5 



Log (Lifetime SFR [M 0 yr 1 ]) 


Figure 4. Recent (t< 200 Myr) SFRs compared with average 
lifetime SFRs for the SHIELD sample. To provide further context, 
the scale on the top x-axis is stellar mass based on the lifetime 
average SFRs and the age of the Universe and assuming a 30% 
recycling fraction. Arrows denote cases where upper limits were 
placed on the SFRs. The solid line represents equivalent recent 
and lifetime SFRs. Overall, the recent SFRs are comparable to, or 
slightly higher than, the lifetime average SFRs. The birthrate pa¬ 
rameters fall in a narrow range with a mean of 1.4 and a dispersion 
of 0.7. The level of recent star-formation activity in each system 
does not appear to be correlated with distance from the nearest 
known galactic neighbors or local environments. 


log (M./MJ 
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Log (Lifetime SFR [M 0 yr 1 ]) 


Figure 5. An expanded view of Figure [4] w ith re sults from 3 
different surveys over-plotted. From fOtt et al 3 mm, typical un¬ 
certainties for the VLA-ANGST recent SFRs are of order ±0.2 
dex; typical unc ertainties for s tellar mass measurements were not 
provided. From [Huang et all Q2012 ). typical uncertainties are of 
order ±0.3 dex for the recent SFRs, and ±0.2 dex for the stel¬ 
lar mass measurements. The shaded region for the SDSS survey 
covers the general range of properties that overlap with the other 
surveys plotted. The recent and lifetime star-formation measure¬ 
ments are consistent with those measured by these three surveys 
in the overlapping mass range. Further, the dTrans are indistin¬ 
guishable from the dlrr galaxies. These low-mass galaxies follow 
the general trends that lower-mass galaxies have lower SFRs. 
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to intrinsic scatter in star-formation activity, we fit a lin¬ 
ear relationship to the recent and lifetime SF Rs using the 
IDL mpfitexy routine (|Williams et al.1 120101) . The mpfi- 
texy routine depends on the mpfit package (jMarkwardtl 
1200911 . We find that less than 35% of the dispersion is 
due to intrinsic scatter while greater than 65% is due to 
the measured uncertainties. Overall the SHIELD sample 
shows recent star-formation activity that is comparable 
to, or slightly higher than, the lifetime average SFR. Fur¬ 
ther, none of the galaxies show star-formation activity 
that is significantly lower than their lifetime average (i.e., 
no quenched star formation), as might be e xpected from 


their relatively low-density environments (|Hodgd 


Einasto et al 


McConnachie 


1974 : Ivan den Berghlll994al lbl: lMated 
2012t IGeha et all 1 201 211. 


1971 


1998 


Figure [5] expands the plot from Figure 0 to include 
results from three other surveys. We add the recent 
SFRs a nd stellar masse s from both the VLA-ANGST 
sample (lOtt et al.l I2012T) and t he sub-sample of ga s-rich 
ALFALFA dwarf galaxies from lHuang et all (l2012f h and 
highlight the gener al range of properties mea sured in the 
SDSS survey from iBrinchmann et all ( 2004 ) that over¬ 


laps with the parameter space in Figure El All stel¬ 
lar masses were normalized to a Salpeter IMF to fa¬ 
cilitate the comparison. Note however, that the stellar 
masses from these surveys were derived using very differ¬ 
ent t echniques including mass-to-ligh t ratios (lOtt et al 


2012), stellar absorption-line indices (iBrinchmann et al 


200411. a nd spectr al energy distribution (SED) fitting 
(iHuang et aljfaofll ). In addition to the ~ 50% uncertain¬ 
ties in stellar mass from these various techniques, there 
are also systematic uncertainties between techniques that 
are not considered. However, these uncertainties are less 
than the broad parameter space probed in FigureEl mak¬ 
ing the comparison between the surveys of value. The 
recent SFRs have also been measured using d i fferent 
techniques. The SFRs from IBrinchmann et ahl (|2004f ) 
were based on emission line model fits and integrated 
Ha e miss ion m easur ements scaled to a SFR. The SFRs 
from lOtt et all ( 2012 ) were derived by scaling integrated 
ultraviolet luminosities. These UV-based SFRs benefit 
from tracing the star-formation activity on comparable 
ti mescale to our CMD-based recent SFRs. The SFRs in 
IHuang et al.l (l2012l l were derived from SED fitting and 
were shown to agree reasonably well with SFRs derived 
by scaling integrated ultraviolet luminosities in systems 
with low extinction. 

There are a number of correlations to note in Fig¬ 
ure El First, while the dTrans (unfilled black points) 
are at the low-mass end of the SHIELD sample, their 
star-formation properties are consistent with the rest of 
the SHIELD sample and are indistinguishable from the 
properties of the dlrrs galaxies measured from the other 
surveys. Second, the trend that gas-rich low-mass galax¬ 
ies generally have higher recent SFRs compared with 
lifetime average SFRs (i.e., b > 1) is more apparent 
across the combined samples. Third, while the properties 
of the SHIELD galaxies are consistent with these other 
studies, the SHIELD galaxies have a narrower range in 
birthrate parameters , parti cularly compared to the AL¬ 
FALFA dwarfs from IHuang et al.l (120121) . One possible 
explanation is tha_t_ the br oader range in the birthrate 
parameter in IHuang et al.l (|2012D is due to higher uncer¬ 
tainties in comparing the stellar masses and SFRs by the 
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Figure 6. The birthrate parameter is pl otted agains t the distance 
to the nearest known neighbor from lMcQuinn et ali >2014 ). While 
the overall range in the b parameter is rather narrow, the galaxies 
with the closest known neighboring system have the lowest values 
of the birthrate parameter, while the galaxies that are located at 
larger distances from their nearest neighbor have higher values of 
the birthrate parameter, albeit with large uncertainties. 

aforementioned different techniques. On the other hand, 
the SHIELD sample is small; future work will expand the 
SHIELD sample to ~30 galaxies and may help determine 
whether the narrow range in recent-to-lifetime SFRs is 
typical in this very low-mass regime (K. B. W. McQuinn 
et al. in preparation). Note that the opposite trend 
(i.e., l ower recent to historical average SFRs) was found 
in iSkillman et al.l (120031) for a sample of gas-rich dwarf 
galaxies in the Sculptor group. However, the SFRs were 
based on integrated Ha emission measurements which 
may not accurately represent the recent star-formation 
activity. SFRs based on Ha emission trace only the 
most recent (~5 Myr) star-formation activity and de¬ 
pend both on constant levels of star formation and on a 
fully-populated upper-end of the IMF. In this low SFR 
regime, not only can the star formation fluctuate on short 
timescales, but the IMF can be stochastically populated 

dBoselli et al.ll2009HLee et al.ll2009l: IGoddard et al.ll201fll: 

iKoda et al.ll2012D Thus, the Ha-based SFRs can be an 
unreliable tracer of star-formation activity, introducing 
scatter in the measurements and potentially biasing SFR 
measurements low in dwarf galaxies. 

5. RECENT SFRS: INFLUENCED BY ENVIRONMENT OR 
REGULATED INTERNALLY? 

One of the outstanding questions concerning the evo¬ 
lution of relatively isolated, low-mass galaxies is the 
comparative importance of environmentally driven star- 
formation activity versus star-formation activity regu¬ 
lated by internal conditions within a system. The former 
is directly correlated with the number of, and proxim¬ 
ity to, neighboring systems. A full investigation of the 
role environment plays requires a complete census of a 
galaxy’s environment including positional and velocity 
information of neighboring systems and SFHs with fine 
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temporal resolution. However, to first order, by com¬ 
paring the recent SFRs of galaxies derived here with 
the proximity of their known neighboring systems de- 
termined from our previous study (iMcQuinn et, al.ll2014L 
see their Figures 5 and 6), we can probe whether there is 
an obvious environmental factor influencing their recent 
evolution. 

In Figure [6] we present a comparison of the b values 
for the sample with the distances to t he nearest known 
neighbors from IMcQuinn et alJ (|2014f ) . All six galaxies 
whose recent star-formation activity is somewhat higher 
than their lifetime average lie in low-density environ¬ 
ments. Specifically, AGC 748778, AGC 174605, and 
AGC 749237 are truly isolated with no known neigh¬ 
bors identified within a 1 Mpc radius. AGC 174585 and 
AGC 731457 are both located ~ 0.9 Mpc from their near¬ 
est known neighbors, namely the 14+ 19 association and 
the low-mass galaxy DDO 83 respectively. AGC 111946 
is located at the end of a linear structure comprised 
of two known ga laxy associations, the NGC 784 and 
NGC 672 groups (iZitrin fc BrosclilfeOQSl IMcQuinn et al.l 
120141) . While AGC 111946 appears to be part of this 
elongated 4 Mpc structure, it is separated by a distance 
of ~ 1 Mpc from the more tightly clustered galaxies in 
the NGC 672 group. Overall, the isolated nature of these 
six galaxies suggests that the recent star-formation ac¬ 
tivity is not driven by gravitational interactions. Note 
however that we cannot rule out the possibility of possi¬ 
ble undetected gas-poor companions outside of the HST 
fields of view. 

Further, the three galaxies that show somewhat lower 
levels of recent star-formation activity compared with 
their lifetime averages (AGC 111164, AGC 111977, and 
AGC 112521) are part of the galaxy associations of 
the NGC 784 and NGC 672 groups mentioned above. 
This linear structure of galaxies stretches ~ 4 Mpc and 
is thought to trace the local matter distribution. In 
contrast to the location of the aforementioned system 
AGC 111946, these three SHIELD systems lie at closer 
distances to the members of these groups. The near¬ 
est neighbor to AGC 111164 is the more massive dwarf 
starburst galaxy NGC 784 at a distance of ~ 0.14 Mpc. 
The SFR in NGC 784 has been noted to be elevated 
over the last few 100 Myr (IMcQuinn et alJfeOlObjh How¬ 
ever, despite being separated by only ~ 0.14 Mpc, the 
recent SFR in AGC 111164 is half its lifetime average. 
AGC 111977 is located ~ 0.65 Mpc from AGC 112521, 
which is also separated by a comparable distance in the 
opposite direction from four members of the NGC 672 
group, namely NGC 672, IC 1727, AGC 111945, and the 
SHIELD galaxy AGC 110482. Even though these galax¬ 
ies reside in a more populated environment, their recent 
star-formation activity is lower than their lifetime values. 
Note that the Hi masses for these galaxies are on the low 
end of our sample (4 x 10 6 M 0 for AGC 111164 and 
7 x 10 6 M 0 for AGC 111977 and AGC 112521), which 
may provide a partial explanation for the low recent star- 
formation activity. 

Finally, we examine the three galaxies whose re¬ 
cent star-formation activity is approximately equiva¬ 
lent to their lifetime averages, namely AGC 110482, 
AGC 182595, and AGC 749241. The first system, 
AGC 1110482, is located in the galaxy association de- 
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Figure 7. Specific SFRs as a function of stellar mass. In contrast 
to the star-forming properties in more massive galaxies, there is a 
high degree of dispersion in these very low-mass galaxies. This dis¬ 
persion likely reflects small fluctuations in star-formation activity 
superposed on already low average SFRs. 

scribed above, > 0.65 Mpc from other members of the 
NGC 672 group. AGC 182595 is truly isolated with 
no known neighbors with 1 Mpc. AGC 749241 lies in 
a more populated region within 9 other galaxies which 
form a linear structure that extends 1.6 Mpc from end 
to end. Given the distribution of galaxies in this re¬ 
gion, it is likely that AGC 794241 has been influenced 
by the external gravitational perturbations from other 
systems. This is supported by the irregular Hi morphol¬ 
ogy whose centroid is not co-spatial with groun d-based 
optic al imaging of the stellar population (I Cannon et alJ 
120111 see their Figure 2). However, the relatively low- 
level of recent star-formation activity in AGC 749241 in¬ 
dicates that any such gravitational interaction has not 
had a dramatic impact on the star forming properties 
over the past 200 Myr. 

In summary, while the range in b parameters is 
fairly small, the somewhat higher levels of recent star- 
formation activity in these low-density galaxy environ¬ 
ments do not correlate with smaller distances to the 
nearest neighboring galaxies, or vice versa, suggest¬ 
ing that the recent activity is governed by internal 
condit ions a nd not e nvironmentally driven. Similarly, 
iHunter fc Elmegreenl (2 0041 ) did not find a correlation 
between star-formation activity and distances to near¬ 
est neighbors in a larger sample of gas-rich dwarf galax¬ 
ies. However, these authors note the possibility that this 
could be a selection effect as there could be undetected 
low-mass neighbors or Hi clouds near the galaxies with 
higher measured recent SFRs. While we cannot rule out 
the possibility of undetected neighboring systems, par¬ 
ticularly if they are gas-poor, there are no systems de¬ 
tected near by in position or velocity in the ALFALFA 
catalog (lHavnes et al.l 120111) . Further, no nearby sys- 
tems have been fo und in the VLA data of the sample 
(j Cannon et al.l[20Tlh . nor in the HST images which cover 




















































Characterizing the Star Formation of the low-mass SHIELD Galaxies 


9 



log (M./M 0 ) 


• SHIELD dlrrs 

2.0 o SHIELD dTrans 



* 


6.5 7.0 7.5 8.0 

log (M,/M 0 ) 


Figure 8. An expanded view of Figure [7] with results from 3 
different surveys over-plotted. Uncertainties for the SSFR mea¬ 
surements are greater than ±0.2 dex for VLA-ANGST and ±0.3 
dex for ALFALFA; uncertainties for the stellar mass were not pro¬ 
vided for the VL A-ANGST measurements and are of order ±0.2 
dex for ALFALFA. The shaded region from IBothwell et al.l (120091) 
covers the range of points included in their study; the extension to 
lower SSFRs includes both measurements and upper limits. Over¬ 
all, for low-mass galaxies none of the samples show a high degree 
of correlation between the SSFRs and stellar mass. 

a field of view extending to ~ 0.5 Mpc for the farther sys¬ 
tems (|McQuinn et alJ[2014D . 

For the galaxies with higher recent SFRs, an alterna¬ 
tive explanation to star formation driven by gravitational 
perturbations may be that these low-mass galaxies have 
experienced a delayed onset to star formation, similar to 
3 low-mass galaxies found in the Local Group (Leo A, 
Leo T, and Aquarius). In a detaile d st udy of the ancient 
star formation history of Leo A, iCole et all (I2Q07T) found 
that over 90% of the stars in the galaxy were formed 
over the mos t rece nt 8 Gyr. In a similar study of Leo T, 
IWeisz et~ahl ( 2012 ) found ~ 50% of the stars in the galaxy 
were formed over the same timescale. The Aquarius dlrr 
also shows delayed onset of star-formation^ with 90% o f 
the stars forming in the last 10 Gyr (iCole et al.l 120141) . 
While a limited comparison based only on three galaxies, 
these examples of delayed or suppressed star formation 
in low-mass galaxies lend support to the idea of galaxy 
downsizing. Higher recent star-formation activity in low- 
mass galaxies may not always be tied to gravitational 
interactions, but may instead be the result of a combi¬ 
nation of internal properties and conditions governed by, 
for example, gas cooling timescales or reionization, and 
stochastic phenomena such as mergers. 

6. COMPARISON OF STAR-FORMATION WITH GLOBAL 
GAS PROPERTIES 

The recent star-formation activity of the sample can 
also be analyzed by normalizing the SFR by the stellar 
mass (M*) of each galaxy. Figure [7] shows the specific 
star-formation rates (SSFR = recent SFR / M*) as a 
function of stellar mass. There is no apparent correlation 
between the SSFR and the stellar mass of the galaxies. 


Figure 9. Fraction of Hi to stellar mass vs. stellar mass for the 
12 SHIELD galaxies. The stellar mass was determined from the 
lifetime average SFRs of the galaxies as described in the text. The 
H i ma ss was calculated from the Hi fluxes reported in i Gannon et all 
mm . The range in gas fractions is consistent with more massive 
gas-rich dlrrs. 

This is in agreement with previous surveys of low-mass 
galaxies which have noted a high degree of scatter in the 
SSFRs as a function of stellar mass in systems with SFRs 
< 10 ~ 3 M 0 yr -1 (e.g. . IBothwell et ahll2009l [Huang et all 

[Ml). 

Figure [5] expands on Figure [7] by addin g results from 
the VLA-ANGST survey ([Ott et al.l I2012T) and a differ¬ 
ent sub-sample of gas-ri ch dwarf galaxies from the AL¬ 
FALFA survey ijHuang et all 12012 ), as well as adding 
th e range o f galax y pro perties measured in the survey 
bv iBothwell et al.l (1200911 . The scat ter is the greatest in 
the study bv iBothwell et all (120091) . particularly below 
M* < 10 7 M 0 , which includes a number of upper limits 
on galaxies with n o detectab l e star formation. Similarly 
to the SFRs fr om iSkillman et all (120031) . the SFRs in 
IBothwell et al.l (120091) were based on integrated Ha emis¬ 
sion whic h ca n be unreli able i n the low-SFR regime. As 
a result, IBothwell et all (I2009H was unable to determine 
whether the scatter in the SSFR-stellar mass relation 
was intrinsic, or d ue to the Ha emissi on not adequately 
tracing the SFR. iHuang et al.l (|2012l) also find a large 
scatter in SSFRs for low-luminosity galaxies, including 
a similar quenching of star formation for some galaxies 
with M* < 10 8 M 0 . However, since th e SFRs in this 
study were based on SED fitting, IHuang et al.l ( 20121 ) 
conclude that the range in SSFR, including the depar¬ 
ture to low SSFRs for some low-mass galaxies, is real. 
In comparison, the SHIELD and VLA-ANGST samples 
show a narrower scatter in the SSFR-stellar mass rela¬ 
tionship, and none of the galaxies in the SHIELD sample 
show quenched star formation (i.e., none have outlying 
low values of SSFR). 

Figure H] presents a comparison of the fraction of Hi 
gas mass to stellar mass (M^z/M*) as a function of 
stellar mass. These gas fractions were calculated us- 
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Figure 10. An expanded view of Figure [9] with resu lts from the 
FIGGS sample over-plotted. From lBegum et alll {2008 1. uncertain¬ 
ties on Mm are °f order 10%, while uncertainties on the stellar 
masses, which are estimated from mass-to-light ratios for the sam¬ 
ple, are not well-quantified. The SHIELD sample overlaps with the 
FIGGS galaxies that have both lower Hi masses and lower stellar 
masses. 

ing the stellar masses derived from the lifetime SFRs 
and t he Hi masses based on Hi flux measurements from 
iCannon et alJ (|201ll ). The Hi masses and gas fractions 
are listed in Table [2] The gas properties show a broad 
distribution as a function of stellar mass with a general 
trend that gas fractions decrease at higher stellar masses. 
The trend and range in gas fraction agrees with those pre¬ 
viously noted for low-mass g alaxies (e.g., Ivan Zee et alJ 
Il997t fSchombert et al .1 [200 if) . Note that while the frac¬ 
tion of Hi mass relative to stellar mass decreases as a 
function of stellar mass, more massive galaxies in our 
sample still have overall higher Hi masses than do their 
less massive counterparts. It is this overall larger reser¬ 
voir of fuel that helps explain the correlation between 
increasing SFRs with increasing stellar mass in gas-rich 
systems seen Figure[H Figure [TU] expands Figure|l]to in¬ 
clude results from the_Faint Ir regula r Galaxies GMRT 
Survey ('FI GGS: IBegum et all I2008D . We use the Hi 
masses from IBegum et al.l (|2008l ) and estimate the stel¬ 
lar masses using the B-band luminosities of the galaxies 
and assuming a mass-to-light ratio of one. The SHIELD 
galaxies overlap with the low-mass end of the FIGGS sys¬ 
tems (i.e., systems with both lower Hi masses and lower 
stellar masses). 

Another way of comparing the SFR and gas content 
is to normalize both parameters by stellar mass. In Fig- 
ureQTj we present a comparison of the SSFR with the gas 
fraction (from above, M hi / M*). Consistent with the 
results presented in Figure [3 the star-formation activity 
for the SHIELD sample shows a high degree of dispersion. 
Figure [12] expands the plot from Figure HD to include 
the results from the VLA-ANGST survey from lOtt et all 
(2012) and a sub-sample of t he ALFALFA survey with 
overlapping mass ranges from iHuang et ah] (120121) . The 
SHIELD sample lies at the low end of the gas ratio ranges 



Figure 11. Specific SFR as a function of the gas ratios (M#j 
/ M*) for the sample. Similar to Figure IT1 there is little correla¬ 
tion between the gas content and star-formation activity once both 
quantities have been normalized by stellar mass. Despite the scat¬ 
ter, no galaxies show significantly lower levels of SSFRs for a given 
gas ratio. 


probed by these other surveys. This may be a selection 
effect as the SHIELD sample was chosen based on low 
Hi masses, and hence may preferentially probe the lower 
range of gas and stellar mass properties. For gas frac¬ 
tions that overlap with the SHIELD sample, the galaxies 
from the other surveys span a larger range in SSFRs. 
This is simply a restatement of our previous finding that 
the SHIELD galaxies have a narrow range in birthrate 
parameter. The broadening in SSFR is in contrast with 
the tight correlation between SSFRs a nd gas f ract ions 
reported for more m assive galaxies (]Brinchmann et alJ 
12004 iLee et afll2007l) and suggests that the stellar mass 
is not the dominant factor in regulating star formation 
in low-mass galaxies. Instead, results in this low-mass 
regime indicate a star-formation process that is depen¬ 
dent upon both the overall gas mass and gas fraction. 
The star-formation variability in gas-rich low-mass sys¬ 
tems is sensitive not only to the overall galaxy mass, but 
also the star formation history of the galaxy due to in¬ 
trinsically smaller gas reservoirs. 

Assuming the star formation in each system persists 
at the recent levels of activity, the gas content in each 
system can be used to estimate how long a system can 
continue forming stars. This gas consumption timescale 
(■ T gas = M on s / SFR where M ffas = 1.33 xMhi to correct 
for helium; lllobertsl 1 1 9631) assumes a 30% recycle frac¬ 
tion of stellar mass (|Kennicutt et all [ 1 9941) and ranges 
from 2 — 34 Gyr; individual values are listed in Table [2] 
Given the low baryonic content of the galaxies, it is un¬ 
likely that their masses have been significantly altered by 
mergers or gas in-fall. Rather, the longer gas consump¬ 
tion timescales for the SHIELD galaxies indicate that 
these low-mass galaxies are inefficient at turning their 
gas content into stars relative to more massive galaxies. 


7. CONCLUSIONS 
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Figure 12. An expanded view of Figure [IT] with results from 2 
different surveys over-plotted. Uncertainties for the SSFR mea¬ 
surements are greater than ±0.2 dex for VLA-ANGST and ±0.3 
dex for ALFALFA; uncertainties for the gas ratio are not quantified 
for VLA-ANGST and of order ±0.1 dex for ALFALFA. There is a 
high degree of dispersion between the SFR and the Hi masses when 
both are normalized by stellar mass. The SHIELD sample overlaps 
with both surveys at the lower end of their gas ratios. This is not 
unexpected as the SHIELD sample was selected, in part, based on 
a criterion of low Hi masses. 

We have used optical imaging of resolved stellar pop¬ 
ulations obtained from the HST to measure the star- 
formation activity of 12 low-mass, gas-rich galaxies in 
the SHIELD program. The average lifetime SFRs range 
from 1(W 3 ' 5 to 10~ 2 2 M 0 yr _1 . The stellar masses mea¬ 
sured from the resolved stellar populations range from 
3 x 10 6 — 7 x 10 7 M 0 , with gas fractions (M hi / M*) be¬ 
tween 0.16 — 1.6. Overall, the galaxies have recent SFRs 
comparable to, or slightly higher than, their lifetime SFR 
averages; none show significantly lower levels of recent 
star-formation activity compared to historical values. 
The ratio of recent to lifetime average SFRs fall in a nar¬ 
row range with a mean value of 1.4 and dispersion of 0.7; 
this is a much smaller range than repo rted in other stud¬ 
ies of gas-rich, low-ma s s galaxies (e.g..lBrinchmann et all 
12004 lQtt et al.l 120121 IHuang et all 120 lT The larger 
scatter found in these other studies may reflect uncer¬ 
tainties in the measurements or in our comparison (see 
Section 4), or may be intrinsic to the galaxies. Future 
work expanding the SHIELD sample to include SFHs of 
30 galaxies will provide a larger, homogeneous sample in 
which to study the birthrate parameter (K. B. W. Mc- 
quinn et al. in preparation). 

We classify two galaxies in the sample, AGC 749241 
and AG C 748 77 8, as dTrans galaxies based on their gas 
content dCannon et al.l 1201 lH an d weak or non-detected 
Ha emission (lHaurberg. Salzer. fc Cannonll20O l. These 
two systems have star-formation and gas properties con¬ 
sistent with dl rrs both in the SHIELD sample and in pre- 


vious studies (Begum et al. 20081: 

Bothwell et al. 20091: 

Ott et al. 2012: Huang et al.l 2012 

). As both galaxies 


have measurable recent SFRs based on their resolved 


stellar populations, the lack of significant Ha emission 
is more likely a result of either stochastic sampling 
of the upper-end of the IMF or a SFR changing on 
short timescales, rather than a fundamental difference 
in galaxy type or properties. 

We used the 3-D dist ribution of galaxies mapped 
by iMcQuinn et al.1 (120141 ) around the SHIELD sample 
to look for correlations between the environments of 
the galaxies and their recent star-formation properties. 
Overall, we find no correlation between the distance to 
the nearest known neighbor and the recent versus life¬ 
time SFR of a system. The more isolated galaxies tend 
to have higher than average SFRs and the galaxies lo¬ 
cated near other systems have lower than average SFRs. 
While we cannot rule out the possibility of an undetected 
system near the galaxies with higher recent SFRs, no 
other systems have been found in the ALFALFA data 
(IHavnes et al.ll2QTTll . these HST i mages (IMcQuinn et al.1 
l2014f) . nor ancillary VLA data dCannon et al.1 1201 If) . 
Thus, the recent star-formation activity in these rela¬ 
tively low-density environments appears to be governed 
by internal evolutionary processes, or the galaxies in 
question may have experienced a delayed onset to star 
formation similar to the g as-rich low-mass g alaxies Leo A 
dCole et al 11200711 . Le o T dWeisz et al .11 2012 b and Aquar¬ 
ius (IQole et all 12011 . 

The star-formation activity in these very low-mass 
galaxies appears to be both inefficient and intrin¬ 
sically variable. In support of an inefficient star- 
formation process, not only are the recent and life¬ 
time SFRs low despite being gas-rich, but the oxy¬ 
gen abundances measured from optical spectroscopy are 
loWj indic ating an i nefficient chemical evolution process 
(lHaurberg. Salzer. fc Cannonll2014l l. These authors have 
shown that 2 of the 12 SHIELD galaxies are classified as 
extremely metal-deficient galaxies (XMD galaxies, e.g., 
iKunth fe Ostlinl l2000h with oxygen abundances below 
12 + log(0/H) < 7.65. In support of a variable, non- 
deterministic process, we find a high degree of dispersion 
between the SSFR, stellar mass, and gas fraction of the 
sample, although we do not report any significantly low 
values of SSFR for a given stellar mass. This idea of vari¬ 
able or “flickering” star-formation activity in low-mass 
galaxies is further supported by the lack of correlation 
between systems with higher recent SFRs and the dis¬ 
tance to the nearest neighboring system. These results 
suggest that the star-formation process is governed by 
the changing internal conditions that are continuously 
shaped by the star-formation process itself. 

Future work on the SHIELD sample will include a spa¬ 
tially resolved study of the Hi mass densities and a com¬ 
parison of the gas properties to the SFRs, infrared emis¬ 
sion from Spitzer warm-mission infrared, and ultraviolet 
emission from GALEX imaging (J. M. Cannon et al. in 
preparation). Additional analysis of the star-formation 
process from the ground-based Ha observations is also 
forthcoming (N. C. Haurberg in preparation). Further, 
our results will be expanded to include the study of a 
larger sample of low-mass galaxies identified from the 
ALFALFA catalog in the SHIELD II survey. 
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